NASA AVSCOM 

Technical Memorandum 105653 f Technical Report 91-C-027 



Optimum Design of a Gearbox 
for Low Vibration 




Katsumi Inoue, Dennis P. Townsend, 
and John J. Coy 
Lewis Research Center 
Cleveland, Ohio 


(NASA-T^-105653) OPTIMUM DESIGN OF A *92-28476 

GEARBOX FOR LOW VIBRATION (NASA) 10 p 


Unci as 

G3/37 0099226 


Prepared for the 

Sixth International Power Transmission and Gearing Conference 
sponsored by the American Society of Mechanical Engineers 
Scottsdale, Arizona, September 13-16, 1992 





SYSTEMS COMMAND 



E-6622 


OPTIMUM DESIGN OF A GEARBOX FOR LOW VIBRATION 

Katsumi Inoue 1 
Tohoku University 
Sendai, Japan 

and 


Dennis P. Townsend and John J. Coy 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

A computer program was developed for designing a low vibra- 
tion gearbox. The code is based on a finite element shell analysis 
method, a modal analysis method, and a structural optimization 
method. In the finite element analysis, a triangular shell element 
with 18 degrees-of-freedom is used. In the optimization method, the 
overall vibration energy of the gearbox is used as the objective 
function and is minimized at the exciting frequency by varying the 
finite element thickness. Modal analysis is used to derive the 
sensitivity of the vibration energy with respect to the design 
variable. The sensitivity is representative of both eigenvalues and 
eigenvectors. The optimum value is computed by the gradient 
projection method and a unidimensional search procedure under the 
constraint condition of constant weight. 

The computer code is applied to a design problem derived from 
an experimental gearbox in use at the NASA Lewis Research 
Center. The top plate and two side plates of the gearbox are rede- 
signed and the contribution of each surface to the total vibration is 
determined. Results show that optimization of the top plate alone is 
effective in reducing total gearbox vibration. 

INTRODUCTION 

The vibration of a power transmission gear system is initiated 
by the gear mesh. The vibration energy is transmitted to the gear- 
box housing through the shafts and bearings, and the gearbox 
housing radiates the structure-borne noise. The vibration and gen- 
erated noise are influenced by manufacturing precision, alignment 
errors of the gears, and the design of the gearbox. These factors 
affect the total quality of the system and are becoming more 
important. 

In order to estimate the vibration of a gear system, it is neces- 
sary to understand the excitation characteristics and to evaluate the 
dynamic behavior of the system. With respect to excitation char- 
acteristics, Kubo, Kiyono, and Fujino (1986) proposed the concept of 
the total vibration exciting force that integrates the effects of 


manufacturing and alignment errors and periodical change of tooth 
mesh stiffness. However, the evaluation of the system’s dynamic 
behavior is related to the modeling of the system. Some approaches 
presented recently include the analysis of multistage gear systems 
coupled with gearbox vibrations (Choy et a!., 1991) and the analysis 
"of vibration transmission in a single-stage gearbox by the building 
block approach (Takatsu et al., 1991). These are not simple simu- 
lations. However, once the vibration of a gearbox is analyzed, the 
sound radiation from the gearbox can be estimated by using the 
results as input data to BEMAP (boundary element method for 
acoustic prediction) (Seybert and Wu, 1989) or a similar boundary- 
element-based, acoustic-intensity analysis. Consequently, the 
vibration and noise of a gearbox can be estimated to some extent. 

Generally, the vibration of the gear system is sensitive to the 
dynamic characteristics of the gearbox. The reduction of noise levels 
is linked to reduction of the gearbox vibration levels. In contrast to 
progress made in the dynamic analysis of gear systems mentioned 
previously, the systematic procedure for reducing gearbox vibration 
levels is not adequate. The technique to avoid a resonance by shift- 
ing the natural frequency is often used to reduce vibration. A mass 
or stiffener is added at an efficient position on the gearbox. Al- 
though this can be an effective method, it is not easy to determine 
the best position on the gearbox or to estimate a suitable natural 
frequency shift. From this viewpoint, the authors have proposed a 
method to design plates (Inoue, Kato, and Ohnuku, 1990) and thin- 
plate structures (Inoue, Townsend, and Coy, 1992). The overall 
vibration energy is adopted as the parameter for estimating the 
quality of a design and is minimized at a given frequency by varying 
the finite element thickness. The optimization problem is solved 
under the constraint of constant weight. In this paper the method is 
applied to the design of a gearbox made from steel plates, and the 
vibration reduction is demonstrated analytically. The exciting forces 
applied at the bearing locations are assumed to be known, although 
influenced by the gearbox’s dynamic characteristics determined in 
the design process. The purpose of this research was to design a 
low-vibration gearbox. 


1 National Research Council - NASA Research Associate. 


OPTIMUM DESIGN PROCEDURE 

Minimization of Overall Vibration Energy 

The theoretical basis of the proposed method is given by Inoue, 
Townsend, and Coy (1992) and is summarized here. The equation 
of motion of a structure is given in the matrix form as follows: 


M{fi(t)}+ [Cj{u(t)} + [K] (u(t)} = {f(t)} 


( 1 ) 


where [M] is the mass matrix, [C] is the damping matrix, [K] is the 
stiffness matrix, {u(t)} is the time- varying displacement vector, and 
(f(t)} is the time-varying exciting force vector. Proportional damp- 
ing is assumed, and the damping matrix is obtained by the linear 
combination of the mass matrix and the stiffness matrix: 

[C] = a[M] + j9|K] ( 2 ) 

where ot and ft are the coefficients to define the proportional 
damping matrix. The modal matrix [4>] is defined by the assembly 
of normalized eigenvectors {^ r }: 


l$] = } , ■ ■ ■ , {<£„} 


(3) 


where n indicates the number of eigenvectors considered. When 
the structure is excited by a harmonic force of angular frequency w, 
the force vector and vibration displacement vector are given by 


fat)} = {F}e^t 


(4) 


where {F} is amplitude of exciting force, {U} is amplitude of dis- 
placement, {£} is the vector defined as {U} = {$}{£}, and j is the 
imaginary unit. Substituting Eq. (4) in Eq. (1), and premultiplying 
it by [$] , the n independent algebraic equations are obtained to 
evaluate the vector {£}. The rth component is given in the fol- 
lowing form by solving the equation 


£r - W T {F} / ((A r - * 2 ) + ju (a + 0X t ) ) 


1, 2, a) 


(5) 


In the expression, X T denotes the rth eigenvalue. Therefore, the 
amplitude {U} for the vibration displacement is obtained by sub- 
stituting {£} into Eq. (4). 

The vibration velocity {u(t)} is given by the similar form as 
shown in Eq. (4). 


{u(t)} = (u}e j “* = m 

The velocity can also be derived directly by differentiating the 
displacement in Eq. (4); therefore, the following equation is ob- 
tained by equating the differentiation to Eq. (6): 

if} = jw {f} 


Therefore, the real and imaginary parts of the rth component ( T 
are obtained as follows: 


frR = (« + ^r)W T fa}/ 


1 (x t - jy + (a + 


( 8 ) 


£rl - - ( A r “ y2 ) |^r} T i F ) ! Hr ( A r ~ + (“ + 0 A r)‘ 

w [w 2 


where R and I are the real and imaginary parts and de- 

notes transpose matrix. The vibration energy T of the structure is 
given by 

t = i (u f m {u} 




(9) 


where (U *} and *} show the conjugate complex of {U} ajid {£}. 
The energy is, therefore, represented by the following expression: 


T = jE (4 + £) 


(JO) 


The sensitivity of the vibration energy with respect to the design 
variable x { (i = 1, . . . , m) is given by 


£T _ A 


£ '» dJ sr + £ « t 

dx I dx: 


( 11 ) 


Since a, ft. w, and {F} in Eq. (§) are independent of the design 
variable, the sensitivity of the vibration energy is represented by a 
function of the sensitivities of the eigenvalue and the eigenvector.. 
The evaluation of the sensitivities of both eigenvalue and eigenvector 
is based on the method proposed by Fox and Kapoor (1968). 

The optimum solution is obtained by minimizing the vibration 
energy under the constraint of constant weight. The element thick- 
ness is adopted as the design variable and limited by specified mini- 
mum and maximum values. The gradient projection method 
(Rosen, 1961) was used for the optimization procedure in this 
research. The value of the energy at every iteration is saved and 
used to determine convergence. The convergence criteria is 1/1000 
relative variation. 


(6) Optimum Des ign Program 

A computer program was developed based on the procedures 
previously described. The program consists of two main pans. The 
first part is the analyzer, which js composed of the finite element 
analysis, the eigenvalue analysis, and the evaluation of the vibration 
energy using modal analysis. The second part is the optimizer, 

(7) which includes the sensitivity analysis and search procedure for the 
optimum value. 


In the finite element analysis, the program is coded to make the 
design program as compact and flexible as possible. If the totcil 
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degrees-of- freedom is constant, using a small number of high- 
precision elements generally leads to good results. However, the 
acceptable degree-of-freedom of the element should be within limits 
in this research because several elements are necessary to model the 
three-dimensional complex surface. Therefore, a triangular shell 
element with 18 degrees-of-freedom is used. The triangular shell 
element is formed by combining the elements for plate bending and 
plane problem analyses given by Zienkiewicz and Cheung (1967). 
The consistent mass matrix used for this analysis is evaluated using 
the numerical integration method. The accuracy of the finite 
element analysis is briefly discussed in the Appendix. 

The program is configured on a Digital Equipment Corporation 
VAX cluster at NASA Lewis Research Center and executed inter- 
actively for small models. For large models the program is run on 
the Cray XMP/28 supercomputer. The program is written in 
FORTRAN 77 with no special software or hardware requirements. 


MODELING AND EIGENVALUE ANALYSIS OF GEARBOX 

Figure 1 shows the gearbox used as a design model. This 
gearbox is used in the NASA Lewis gear-noise test facility (Oswald 
et al., 1991) and is made of steel plates. The origin of a Cartesian 
coordinate frame is placed at the center of the bottom plate, and the 
positive Z axis is defined upward and normal to the bottom plate. 
The Y axis is parallel to the direction of the shafts. The overall 
dimensions of the gearbox are approximately 324 by 248 by 270 mm 
(dimensions in X, Y, and Z directions). The four side plates are 
welded together and to the bottom plate. The top plate is bolted to 
the side plates. To simplify the formulation, the gearbox is modeled 
as a simple box by using 256 triangular shell elements with all 
surfaces continuously connected. The stiffeners and the four circular 
holes for the bearings are neglected in the model as is the supportive 
effect of the shafts. These simplifications may be reexamined and 
more elements used for a precise analysis, while still conserving the 
basic concept of the proposed method. The number of elements is 
fixed in this research only to limit computer time. The four corners 
of the bottom are fixed. Figure 2 shows the mesh pattern. 




Figure 2. — Mesh pattern used for the optimum design of NASA 
gearbox. 


The six surfaces are identified by the numbers 1 to 6 as follows: 
surface 1, the top; surfaces 2 and 4, the plates at Y — 124 and 
-124; surfaces 3 and 5, the plates at X= 162 and —162; and sur- 
face 6, the bottom. The thickness of the portions for the bearing 
housings at -121 < X < 82 and 66.5 < Z < 193.5, is 22.4 mm; the 
other elements are 6.3 mm thick. The constants used in this re- 
search are modulus of elasticity ~ 206 GPa, Poisson's ratio = 0.3, 
mass density — 8000 kg/m 3 , coefficients for composing proportional 
damping matrix a = 1.0 sec* 1 , and = 5.0x10 sec. 

The eigenvalue problem is solved first to obtain the dynamic 
characteristic of the gearbox. The lower 10 natural frequencies arid 
mode shapes are summarized in Table X. The mode shape of each 
surface is represented as, for example, (m x , m^), which indicates that 
n^— 1 and m^— 1 nodal lines are observed in X and Y directions. 
The calculated fundamental natural frequency of 475 Hz is approxi- 
mately 4.4 percent less than the result of measurement in the report 
by Lim, Singh, and Zakrajsek (1989). Some other measurement 
data are shown in this report, but there are no modes that coincide 
with the calculated results. This may be caused by neglecting the 
effects of stiffeners and shafts. 


VIBRATION OF GEARBOX RESULTING FROM UNIT 
HARMONIC FORCES 

Unit harmonic forces IN are applied to the Four positions of the 
bearings in the direction of the line of action, and the total vibration 
energy is calculated for the exciting frequency, which is given at 
every 10 Hz between 500 to 1500 Hz. The energy is evaluated on 
the basis of the lower 30 modes. Since the natural frequency of the 
30th mode is 2448 Hz, which is higher than the considered frequency 
range, 30 modes is enough for this calculation. The vibration energy 
is shown by the broken line in Fig. 5. Ten peaks are observed in the 
range at 540, 580, 610, 770, 940, 1240, 1280, 1320, 1400, and 
1460 Hz, respectively. Since the energy is calculated at every 10 Hz, 
the peaks do not coincide exactly with the natural frequencies, some 
of which are shown in Table I. However, they indicate the 
resonance approximately . 
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TABLE 1. — LOWER 10 NATURAL FREQUENCIES AND MODE SHAPES OF GEARBOX 


Natural 

frequency 

number, 

i 

Fre- 

quency 

(U. 

lie 

Surface number, mode shape 


2(cn A ,m t ) 

3(m y ,in m ) 

4(m xl m,) 

5{m ,mj 

6(m x ,m r ) 

, 

475 

*0,0 

*0.0 

‘(M) 

*0,1) 

‘(1,0 

*(1.0 

2 

540 

(1,2) 

*0,0 

(2,0 

‘0.0 

(2,1) 

(1,2) 

3 

543 

*0.0 

0.2) 

,(..*) 

(1,2) 

(1,2) 

‘(1,0 

4 

582 

‘0.0 

‘0.0 

b 0,0 

‘{1.0 

*0-0 

*0,0 

S 

607 

(2.0 

(2,1) 

‘0,0 

(2,0 

b (l,l) 

(2.0 

6 

768 

*0.0 

*(1,0 

‘0,0 

*0.0 

•0.0 

*0,0 

7 

942 

(2.1) 

C {1,2) 

(1,2) 

‘(1,2) 

(1,2) 

(2.1) 

8 

991 

(2,0 

(2.0 

‘0.0 

(2.0 

*(1,1) 

(2.0 

9 

1041 

(1.2) 

‘(1,2) 

(2,0 

‘0,2) 

(2,0 

0,2) 

10 

noi 

(1.2) 

‘0*2) 

(2.0 

*0»2) 

i (2,0 

(1.2) 


'Convex outward. 
^Concave outward. 
'Oblique nodal line. 


Generally, the product of an exciting force and the correspond- 
ing displacement gives the work done by the force. Since the work 
done represents the energy input into the vibrating structure, it 
must have a strong correlation with the vibration energy defined in 
this research. This correlation is demonstrated approximately for 
the aforementioned frequencies of the gearbox in the following 
manner. 


The mode shape obtained by the eigenvalue analysis is not the 
actual displacement, but the relative displacement at every node. 
However, all the eigenvectors calculated in this research are normal- 
ized with respect to the mass matrix. Therefore, if every eigenvector 
component is listed at a certain node, it might represent the relative 
linear and rotational components of the displacement at the node. 
The components at the position of the bearings are selected from the 
results of the calculation and normalized with respect to the maxi- 
mum value to obtain the displacements. Figure 3 illustrates the dis- 
placements at bearing I. The abscissa indicates the mode order 
from 1st to 30th. For example (Figs. 3(a) and (c)), the maximum 
displacements in the X and Z directions originate in the 27th and 
28th mode shape, respectively. 

The product of the X and Z components of the applied unit 
load and the corresponding displacements at the bearing positions is 
defined here as the work done by the load, 



(a) Displacement along X axis and 
rotation about Z axis. 


Y-displacement 
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(b) Displacement along Y axis. 


U ~ |uj - u 2 | sin a w + Jwj - w 2 | cos (12) 

where u and w are the displacements in X and Z directions; 
a w is the operating pressure angle; subscripts 1 and 2 indicate the 
position of bearings 1 and 2, respectively. The displacements 
and Wj are taken from the results at the aforementioned 10 fre- 
quencies: the 2nd, 4th, 5th, 6th, 7th, 11th, 12th, 13th, 15th, and 
17th modes in Fig. 3. A similar method is used to obtain u 2 and 
w 2 , and the work done is calculated from Eq. (12). The work done 
is not the actual, but the relative work. Figure 4 shows the com- 
parison of the work and the calculated vibration energy of the gear- 
box. The energy is approximately proportional to the work. This is 
as expected and suggests that the relative vibration magnitude can 
be estimated approximately from the work. 

OPTIMUM DESIGN OF GEARBOX 
Design of the Top Plate 



The vibration energy of the gearbox caused by unit exciting 
force IN with initial thickness is shown by the broken line in Fig. 5. 
The abscissa indicates the exciting frequency, and the peaks repre- 
sent the resonance. In the first example, the top plate of the gear- 
box is designed through the energy minimization at 1260 Hz. This 
frequency coincides with the gear mesh frequency for a 28-tooth gear 


(c) Displacement along Z axis and 
rotation about X axis. 

Figure 3. — Relative displacements at the 
positions of bearing 1 (normalized by the 
maximum displacement). 
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Work done by excitation 

igure 4. — Relation between vibration energy of gearbox and 
work done by exciting force. (Work done is calculated from 
the normalized displacements.) • 



500 1000 1500 

Frequency, Hz 

Figure 5.— Frequency response of gearbox due to unit exciting 
force. (Top plate is optimized.) 


operating at 2700 rpm. The initial vibration energy at this fre- 
quency is about 0.976xl0‘ 8 J. The limits of thickness are specified 
as 2 mm < t < 12 mm. 

The energy converges to a value of about 19 percent of the 
initial value after six iterations (Fig. 6). The frequency response 
after optimization is shown by the solid line in Fig. 5. The 
resonance peak at 1280 Hz shifts to approximately 1320 Hz. The 
response variation may look slight in the figure, but the vibration 
energy is reduced. This optimum design solution required about 
875 sec per iteration on the Cray XMP/28. 



Iteration number 

Figure 6.— Convergence of vibration 
energy. (Top plate is optimized.) 


Figure 7 illustrates the optimum shape of the top plate. The 
thickness increases above bearing 1 and near the edge that connects 
to surfaces 2 and 4. The displacement along the Y axis at bearing 
1 is large in the 12th mode, 1277 Hz (Fig. 3(b)). The increase in 
thickness causes an increase of rotational rigidity about the X axis 
to control the displacement. Figure 8 shows the energy share of 
each surface before and after the optimization. The energy reduction 
in the top plate is comparatively small. However, the variation in 
thickness of the top plate causes a vibration energy reduction in the 
other plates. 


Design of the Top Plate and Two Side Plates 


In the second example, two side plates 3 and 5 as well as the 
top plate are optimally designed^ The side plates are selected 
because they do not support the shafts. The conditions of excitation 
and constraint are the same as in the first example. Figure 9 shows 
the frequency response. The energy is considerably reduced between 
1180 to 1370 Hz, and the vibration energy reduction is larger than in 
the first example. The energy is reduced sharply at the first itera- 
tion and converges on the value of approximately 3 percent of the 
initial value (Fig. 10). 



Figure 7.— Optimum shape of gearbox top plate. 
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Figure 10. — Convergence of vibration 
energy, (Top plate and two sides 
are optimized.) 



Figure 9.— Frequency response of gearbox due to unit exciting 
force. (Top plate and two sides are optimized.) 


The optimum shape of the top plate and two side plates are 
illustrated in Fig. 11. The thickness distribution of the top plate is 
similar to the result in the first example. However > the weight of 
the top plate is reduced approximately 15 percent by the optimiza- 
tion as compared to the initial weight. The optimum shapes of the 
two side plates are similar to each other. The thickness increases at 
the bottom and sides of the plates. The weight of plate 5 is also 
reduced approximately 7 percent from the initial weight. However, 
the weight of plate 3 increases approximately 23 percent because of 
the constraint condition of constant weight. Figure 12 shows the 
energy share of each surface. The energy is reduced in every surface, 
and the reduction in plate 3 is substantial. This may be partly due 
to the increase in weight. 

Verification and Revision of the Proposed Method 

An experimental verification may be required for this kind of 
design problem for two main reasons: to verify the proposed method 
and to estimate the error of the result. Method verification can be 
done using computer simulation instead of experimentation. After 
obtaining the optimum shape of a beam, the authors showed that 
the displacement amplitude was reduced considerably f ollo wing opti- 
mization by solving the forced vibration problem of the beam (Inoue, 
Townsend, and Coy, 1992). This validates the propose d m e tho d it- 
self. Estimating the error of the result depends on the assumptions 
and simplifications of modeling discussed in the section on modeling 
and eigenvalue analysis. A comparison with the measured frequency 
is presented also in this section. 

A revision of the method could be completed later and applied 
to the design of low vibration and lightweight gearboxes by elimina- 
ting the constraint of constant weight and introducing appropriate 
side constraints. 
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Energy, J 


CONCLUSION 




Figure 11. — Optimum shape of the top plate and two side plates 
of gearbox. 



Gearbox face 


Figure 12. — Share of vibration energy. (Top plate and two side 
plates are optimized.) 


The method proposed to design a low vibration gearbox is based 
on a finite element shell analysis, modal analysis, and a structural 
optimization technique. The overall vibration energy is adopted as 
the objective function and minimized by varying the element thick- 
ness under the constraint of constant weight. A computer program 
was developed for this purpose. 

The computer program was applied to the design problem of a 
test stand gearbox used at NASA Lewis. The gearbox was excited 
by unit harmonic forces applied at the bearing positions in the 
direction of the line of action. The work done was estimated by the 
product of the bearing forces and the relative displacements obtained 
from the eigenvalue analysis. The work was approximately propor- 
tional to the total vibration energy of the gearbox. This suggests 
that the gearbox vibration magnitude can be estimated from the 
relative work. 

The top plate and the two side plates of the gearbox were opti- 
mized. The vibration energy was reduced at the excitation fre- 
quency. The vibration reduction was obtained by variation of the 
thickness distribution and the consequent shift of weight from the 
top plate to the side plate. The contribution of each surface on the 
total vibration was determined. Optimization of the top plate only 
had considerable effect in reducing the vibration. 



Figure 13.— Error of eigenvalue analysis of simply supported 
rectangular plate. 
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APPENDIX — ACCURACY OF THE FINITE ELEMENT 
ANALYSIS 


The eigenvalue problem of a simply supported rectangular plate 
is solved by the developed program. The plate is divided into 
4x4x2, 6x6x2, and 8x8x2 elements, respectively. The mesh pat- 
terns are similar to the pattern of the top surface of the gearbox 
shown in Fig. 2. Most of the computed natural frequencies are 
smaller than the analytical values, and the error depends on the 
mode shape or the mode order. The squared error arithmetic mean 
is evaluated and its square root is illustrated in Fig. 13. In the 
lower 30 gearbox modes, which are used for the modal analysis, each 
surface has 1 and 3 nodal lines at most. This mode corresponds to 
the 10th mode of the plate shown in Fig. 13, suggesting that the FE 
mesh in each surface of the gearbox may be used for the calculation 
of the mode shape with an error of approximately 10 percent. This 
does not directly show the error of the gearbox vibration because of 
its complexity. 
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13. ABSTRACT (Maximum 200 words) 

A computer program was developed for designing a low vibration gearbox. The code is based on a finite element shell 
analysis, a modal analysis, and a structural optimization method. In the finite element analysis, a triangular shell element 
with 18 degrees-of-freedom is used. In the optimization method, the overall vibration energy of the gearbox is used as 
the objective function and is minimized at the exciting frequency by varying the finite element thickness. Modal analysis 
is used to derive the sensitivity of the vibration energy with respect to the design variable. The sensitivity is representative 
of both eigenvalues and eigenvectors. The optimum value is computed by the gradient projection method and a 
unidimensional search procedure under the constraint condition of constant weight. The computer code is applied to a 
design problem derived from an experimental gearbox in use at the NASA Lewis Research Center. The top plate and two 
side plates of the gearbox are redesigned and the contribution of each surface to the total vibration is determined. Results 
show that optimization of the top plate alone is effective in reducing total gearbox vibration. 
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